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Abstract. Stable laws can be tempered by modifying the Levy measure to cool 
the probability of large jumps. Tempered stable laws retain their signature power 
law behavior at infinity, and infinite divisibility. This paper develops random walk 
models that converge to a tempered stable law under a triangular array scheme. 
Since tempered stable laws and processes are useful in statistical physics, these 
random walk models can provide a basic physical model for the underlying physical 
phenomena. 
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1. Introduction 

Tempered stable laws were i ntroduced in p hysics as a model for turbulent velocity 

so been used in finance 



CarretaL 



fluctuations ( Koponenl (j!995h : iNovikovl (119941) ). They have a 



(120021 120031 )) and hydrology (IMeerschaert et al. 



transient anomalous diffusion (JBaeumer and Meerschaert 



(120081 )) as a model of 



(120101 )) . The general clas s 



of tempered stable distributions for random vectors was developed by lRosihskil (120071 ). 
In short, the Levy measure of a stable law is modified in the tail to reduce the 
probability of large jumps. Often this is done in such a way t hat all moments exis t, 
but tempering by a power law of higher order is also useful (jSokolov et al.l ( 120041 ) ). 
This paper develops random walk models that converge to a tempered stable law. 
Starting with a random walk in the domain of attraction of a stable law, the basic 
idea is to modify the tails of the jumps to mimic the tempering function of the limit. A 
triangular array scheme is essential, since the limit is no longer stable. The results of 
this paper are intended to form a useful random walk model for natural processes that 
are well described by a tempered stable. The main result of this paper is Theorem 
14.31 which shows that the weak limit of the row sum of that triangular array is a 
tempered stable distribution. In Theorem 14.81 we show that the random walk process 
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2 ARIJIT CHAKRABARTY AND MARK M. MEERSCHAERT 

converges to the Levy process generated by the tempered stable distribution in the 
sense of finite-dimensional distributions. 

Section [2] gives a brief background of stable distributions and their domains of 
attraction. In Section [3l we define tempered stable distributions and the triangular 
array model. In Section HI we state and prove the results regarding the convergence 
of the random walk to the tempered stable distribution. 



2. Stable limits for random walks 

Recall that a random vector X on R d is infinitely divisible if and only if its char- 
acteristic function _E[e^ A ' x ^] = e^ x \ where 

(2.1) V(A) =i(a,\)-l(\,Q\)+ [ (e t{x ' x) - 1 -i(A,x)l(||x|| < l)W(dz), 

where a G M. d , Q is a nonnegative definite symmetric d X d matrix with entries in R, 
and M is a cr-finite Borel measure on M. d \ {0} such that J" ,„min{l, ||x|| 2 }M(<ix) < 
oo. The triple [a, Q, M] is ca lled the Levy representation, and it is unique 



( IMeerschaert and Schefflerl . I2001L Theorem 3.1.11). The measure M is known as 
the Levy measure of X. 

A M. d valued random vector X is said to be stable if and only if for all n > 1, there 
exist b n > and a n G R d so that Xi + . . . + X n = b n X + a n , where X\, X 2 , . . . are 
i.i.d. copies of X. Clearly, a stable random vector is infinitely divisible. It is well 
known that given a stable random vector, either it is Gaussian, or the Gaussian part 
is completely absent, i.e., in the Levy representation ( 12. ip . Q = 0. In this paper, 
"stable random vector" will refer to the latter case, i.e., non-Gaussian stable random 
vectors. For such a random vector X, P(||X|| > •) regularly varies with index —a 
for some < a < 2. Sometimes, X is also referred to as an a-stable random vector. 
The Levy representation of the random vector X is [a, 0, M ] for some a G R d where 
M (dr,ds) = r~~ a ~ 1 dra(ds), and a is a finite non- zero Borel measure on the unit 



sphere S 1 = {x G R : ||x|| = 1}, see for example (IMeerschaert and Scheffler 



2001 



Theorem 7.3.16). The measure a is known as the spectr al measure of X. For more 

detail s on stable distributions, the reader is referred to 

il994T ). 



Samorodnitsky and Taqqu 
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The domain of attraction of an a-stable random vector X consists of M. d valued 
random vectors H such that there exist a sequence of positive numbers (b n ) and a 
sequence (a n ) in M. d satisfying 

(2.2) b- 1 (H 1 + --- + H n ) -a n ^X 

as n — y oo, where =>■ denotes weak convergence and Hi,H2, ■ ■ ■ are i.i.d. copies of 
H. A necessary and sufficient condition for (12. 21) is that V(r) 
regularly with index —a, and 

P(\\H\\:> r.^ G U) ai^D) 



P{\\H\\ > r) varies 



(2.3) 



PiSrr eD\\\H\\ >r) 

\ \\-ti I " M / 



> r >pir e£ 



— y 



V(r) ' a(S d - 1 ) 

as r — > oo for all Borel subsets D of S^ -1 such that cx(dD) = 0, see for example 
(IMeerschaert and Schefflerl . l200ll . Theorem 7.3.16). When a > 1, m — E(H) exists, 
and we can center to zero expectation in ( 12.21) by setting a n = nh~ l m. Then the limit 
X also has zero mean, and its log-characteristic function 

(2.4) V(A) = / U {x ' x) - l-i{X,x))M{dx). 

When < a < 1, no centering is required: Set a, 
stable law with log-characteristic function 

3 J(A,x> 



in (12.21) ; then X is a centered 



(2.5) 



VW 



x^O 



l)M(dx). 



2001 



Theorem 8.2.16). 



See for example ( IMeerschaert and Schefflerl . 

Suppose that X is a stable random vector. Let {X(t)} denote the Levy process 
associated with X, so that X(0) = almost surely, X(t) has stationary, independent 
increments, and X(l) = X in distribution. Suppose that ( 12. 2 p holds with a n = 0, let 
6(c) = 6|- c ], and S(c) = Yljli-Hj f° r c — 0- Then, as c — y oo, we also have process 

{ X(t)}/ >n in the sense of finite dimensional distri- 



2001 



Example 11.2.18) as well as convergence in 



converge nce |6(c)~ 1 >5'(ct)|^>n 
butions ( IMeerschaert and Scheffler 
the Skorokhod space D (\Q, oo), R rf ) of right continuous f unctions with left hand lim- 
its, in the J\ topology (IMeerschaert and Schefflerl . I2004J . Theorem 4.1). The random 
vectors X(t) have smooth density functions P(x,t) that s olve a fractional dif f erenti al 



equation -^P(x,t) = DV"F(i,t) for anomalous diffusion (IMeerschaert et al.l (I1999I )). 
The fractional derivative operator V^ f /(x) is defined, for suitable functions f(x) 
with Fourier transform /(A) = f e l ^ x ' x 'f(x) dx, as the inverse Fourier transform of 
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L |i =1 (— «(A, s)) a a(ds)f(\), and V > is a positive constant that depends on choice 
of the norming sequence b(c). The random walk S n provides a physical model for par- 
ticle jumps, whose ensemble behavior is approximated by the stable density functions 
P(x, t). For example, the random walk can be simulated to solv e the fractional diff u- 



sion equation, a numerical method known as particle tracking (iZhang et al 



(H)). 



The purpose of this paper is to develop analogous random walk models for tempered 
stables. 



3. Tempered random walks 

This section develops random walk models th a t con verge to a tempered stable, 
using a triangular array scheme. As in iRosihskil ( 120071 ). we define a rf-dimensional 
proper tempered a-stable random vector to be an infinitely divisible random vector 
with Levy representation [a, 0, M] with 



(3.1) 



M(dr, ds) = r a q(r } s)dra(ds 



for r > and s G S d ~ x . Here a G (0,2), a is a finite Borel measure on the unit 
sphere S' d_1 = {x G M. d : ||x|| = 1}, and q : (0, oo) x S^ -1 is a Borel measur- 
able function such that for all s G S^ -1 , q(-,s) is non- increasing, q(0+, s) = a and 
q(oo, s) = 0. We also assume that q is c ontinu ous in the second variable, i.e., q(r, ■) 



is continuous for all r > 0. In IRosihskil ( J2007I ). the assumption is that q(0+, s) =- 1. 
However , a sim ple reparametrization yields q(0+, s) = a. It is also assumed in 
Rosihskil ( 120071 ) that q{-,s) is completely monotone, but we do not need that as- 



sumption in this paper. Note that tempered stable random vectors are full dimen- 
sional, since the Levy measure is not concentrated on any lower dimensional subspace 



( IMeerschaert and Schefflerl . |2001| . Proposition 3.1.20). Let if be a random vector in 
the domain of attraction of X such that (without loss of generality) P(H = 0) = 0. 
We will define a random walk that approximates the tempered stable using a condi- 



tional tempering of the jumps. Define a function it : (0, oo) x S' 



d-l 



— )• 



by 



(3.2) 



7T{U, S 



U ' 



r a q(r,s)dr. 
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From the fact that the function q is bounded above by a, it is immediate that the 
integral on the right hand side is finite. Clearly, 

dn{u, s) 



au a 1 r a 1 q(r,s)dr — u 1 q(u,s) 

J u 
/•oo 

a/" 1 / r" a ~ l {q(r, s) - q(u, s)} dr < 0, 



/•oo 
.a— 1 / _— Q— 1 



u 



the inequality following from the fact that q(-, s) is non-increasing. Thus ir(-, s) is also 
non-increasing. A simple application of the L'Hopital's rule yields that vr(0+, s) = 1 
and 7r(oo, s) = 0. 

Define a family of probability measures on (0, oo) by Ii(du, s) = —-^ir(u,s)du. 
Since q is a measurable function, for every Borel set A C (0, oo), 11(74, •) is a mea- 
surable function from 5' d_1 to R. Hence there exists a random variable T > whose 
conditional distribution given H is II (•, H/\\H\\). 

Now we construct the tempered random walk. For t > 0, define 

(3.3) H^^-MHWMT). 

\\ti || 

Let {(Hj : Tj) : j > 1} be i.i.d. copies of (H,T). Suppose v n — > oo is a sequence of 
positive numbers. Define Y n j := (||-Hj|| A v n Tj) Hj/\\Hj\\ for n, j > 1, and let 

(3.4) Sn (k):=Y nl + --- + Y nk . 

The next section shows that, for suitably chosen truncation thresholds v n , the random 
walk (j3.4p is asymptotically tempered stable. 

4. Limits of tempered random walks 

The results in this section show that for suitable truncation thresholds v n and 
centering, the random walk (13. 4 j) converges to a tempered stable. The form of the 
centering is then related to the case with no tempering. We begin with a few technical 
results. Recall that for sigma-finite Borel measures \i n ,\i on T = ~R d \ {0}, \x n — > \i 
(vague convergence) means that fi n (D) — > fi(D) for Borel sets D G T that are bounded 
away from zero, for which n(dD) = 0. 

Lemma 4.1. Suppose that H l is defined by (13. 3p . Then 

P(||/f || > *) fffS- 1 - 1 ) v ' 



6 ARIJIT CHAKRABARTY AND MARK M. MEERSCHAERT 

on T as t —y oo ; where M is as in ( 13.11) . 



For the proof, we shall need t he following; result from wea 



is similar to Theorem 1.3.4 in 



van der Vaart and Wellner 



£ conve rgence. This result 
(Il996f ). which states the 



corresponding result for probability measures. Although the result is well known, we 
include the proof here for completeness, since we could not locate a suitable reference. 

Lemma 4.2. Suppose (/i„) is a sequence of measures on some metric space S 
equipped with the Borel sigma-field, converging weakly to some finite measure fi. 
Then, for all bounded non-negative upper semicontinuous functions f , we have 
limsup^^ / fdfin < J fdfi. 

Proof. Since / is bounded and non-negative, we can assume without loss of generality 



that < / < 1. For k > 1, denote f k := J Ya=i 1 {f~ X [ir, l))- It is easy to 
that / < fk < / + h- Thus, for fixed k > 1, 



sec 



n— >oo 
k 



lim sup / fd/j, n < lim sup / f k d/j, n < - ^ lim sup fi n ( f 1 



j=i 



i- 1 



<-lE,(ri^))<-I 



! ^,i)) <//*+>), 



by the Portmanteau Theorem (Theorem 2.1. iBillingsleyi (119681 )) and the observation 
that / _1 [M^, l) is a closed set because / is upper semicontinuous. Since \x is a finite 
measure, this completes the proof. □ 



Proof of Lemma \4-1\ Note that if q(r,-) is continuous for all r > 0, then the same 
is true for tt(u, •) for all u > 0. To see this, suppose that the former holds. Fix 
u > and s n ,s E S d ~ l such that s n — > s. By the assumption on q and the dominated 
convergence theorem, it follows that lim^^oo J r~ a-1 g(r, s n )dr = J r _a_1 g(r, s)dr . 
Then ii(u, s n ) —¥ n(u, s), so tt(u, •) is continuous for every u > 0. A similar argument 



shows that c,7r ^'' j is continuous for every u > 0. 



(4.1 



dn(u,-) 
du 

For the proof, we shall use the fact that as t — y oo 

H 



P(\\H\\ >ty 1 P 



t 1 \\H\\ e dr, -—— e ds 



H\ 



— y 



a 



-a-l 



a(S d 



-V 



dra(ds) 
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which is a restatement of Theorem 8.2.18 in 



Meerschaert and Schemer! (120011 ). It 



suffices to show that for every closed set Act \ {0}, 

(4.2) limsupP^P* G A)P{\\H\\ > t)- 1 < a(S d - 1 )- 1 M(A) , 

t— >oo 

and that for every e > 0, 

(4.3) lim P{t- l H l G B C £ )P(\\H\\ > t)' 1 = a(5^ 1 )- 1 M(5 £ c ) , 

t— >oo 

where B r is the closed ball of radius r centered at origin: B r := {x G M d : ||x|| < r} . 
Fix a closed set A C R d \ {0} and note that 

/*0O P [*CG TT 

P(r 1 H t eA)= / / l A ((rAu)s)U(du,s)P(r 1 \\H\\edr,-—eds) 

JO J 3 d - 1 Jo II-" II 

Fix sequences r n and s n so that r n > 0, r n — > r > 0, s n G S^ 1 and s n —¥ s. Then 

/■oo 

lim sup / 1a ((r n Au)s n )Tl(du, s n ) 

n— >oo Jo 

"^ ' dn(u,s r 



< 1- / liminf l j4 c((r„ Au)s„) — du 

Jo n -*°° \ ou J 

f°° ( dn(u s)\ f°° 

= 1-1 )-l—L )]im.mfl A c((r n A u)s n )du< / l A ((r A u)s)IL(du, s) 

Jo V * / n ^°° Jo 

by Fatou's Lemma, continuity of J^7r(w, •), and the fact that A c is open. Then (r, s) i— >■ 
J" 1a((?" A u)s)H(du, s) is upper semicontinuous. From ( 14. ip . it follows that for all 
e > 0, the restriction of P(||iP|| > t) _1 P [t _1 ||iJ|| G dr, H/\\H\\ G ds] to P £ converges 
weakly to that of ar~ a ~ 1 dra(ds)/a(S d ~ 1 ). Thus, by Lemma [4.21 and the fact that A 
is bounded away from zero, it follows that 

*Z!? %mt?) £ [ L [ 1Mrh u)s)Jl(du ' s ^f a - ldr °^ 
= 1°° L {[ lA(vs)n{du - s) } Tuh-f^^ 

+ [ /{/ lA(rs)n(du,s)\ ° r-'^dMds) =: h+h. 

Note that Lemma [4.21 applies since U(du, s) is a probability measure for each s. 
A change of the order of integration yields that 

h = ^i) / L-, lA{us)u ~ a (-^r ) aWdu ■ 
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It is immediate that li = cr(S' d_1 )~ 1 J °° f sd -i lA(rs)a7r(r, s)r~ a ~ l dra(ds) . Thus, 
P^H 1 E A) 

lim SUP —77! — J, — 

(44) <->J P(\\H\\>t) 



- ^S^)J J sd _Urs)^(r,s)-r^-jr— dra(ds) . 

From (|3D, it follows that «7r(r, s) — r^7r(r, s) = g(r, s) . Plugging this in (j4.4p yields 
that 

Pif- 1 !!* e A) 1 [°° f _ a-1 _ M (A) 



hmsup ^7T ^ ^ed-n / / U{rs)q{r,s)r dra{ds) - , 

thus showing (14.21) . For (14.31) . note that as t — > oo, 
Pit' 1 !! 1 e B c £ ) = f f 7r(e, s)P (V 1 ^!! 6 dr, t-^ G ds ) 

A ^-1 V II-" II / 

~ P(\\H\\ >t) [^ f 7l(e, s)-^r— 1 draids) 

by (14.1 p . the fact that 7r(e, •) is continuous, and the definition of n. D 

The following theorem is the main result of this paper. It shows that the tempered 
random walk (j3.4p converges weakly to a tempered stable, for suitable tempering 
constants v n and suitable centering vectors a n . 

Theorem 4.3. Forn > 1 letb n := inf {x : P(\\H\\ > x) < n^ 1 } . If the sequence (v n ) 
satisfies 

(4.5) limt;- 1 6 n = «T(S , ^ 1 ) 1 / a , 

then, 

(4.6) v~ x S n {n) -a n ^> p 

where p is an infinitely divisible probability measure on M. d with !evy representation 
[0, 0,M], and (a n ) is defined by 



(4.7) a n := n 



f xP{v- x Y nl e dx) 

J{\\x\\<\} 
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Proof. Note that nP(\\H\\ > b n ) -> 1 (IResnickl . 120071 p. 24). Since P{\\H\\ > •) varies 
regularly with index -a, (1451) implies P(||-ff|| > v n )/P(\\H\\ > b n ) -)• a(S' d_1 ) . Then 



<t(5' 



(4.8) lim nP( \\ H || > v n 

n— >oo 

An appeal to Lemma 14.11 shows that 

(4.9) nP{v~ n x Y nX G •) A M(-), 



-i> 



and then it suffices to check (iMeerschaert and Scheffler 



2001 



(4.10) 



lim lim sup nv n 2 E ||F„il(||F„i|| < v n l 



Theorem 3.2.2): 
0. 



For this, note that ||y n i|| 2 l(||y nl || < vj) < ||Pi|| 2 l(||i/i|| < vj) + u 2 T 2 l«Ti < 
V, ll^xll > vj) < WHtfliW^W < vj) + v^liWH.W > vj). Since P(\\H\\>-) is 



2007 



regularly varying with index —a and a < 2, by Karamata's theorem (jResnick 
Theorem 2.1) it follows that E [\\H\\ 2 1(\\H\\ < vj)} ~ K5) 2 P(||P|| > u n 6)a/(2 - 
a) ~ 5 2 ~ a t> 2 P(||P|| > f„)a/(2 — a) as n — )• oo. Using the regular variation of 
P(||iP|| > •) once again, it is immediate that E [v^8 2 l(\\H\\ > v n 5)} = i; 2 <5 2 P(||P|| > 
v n 5) ~ 5 2 ~ a v 2 l P(\\H\\ > v n ) . To complete the proof, use (14.81) to obtain C < oo such 
that for all 5 > 0, 



(4.11) 



limsupm;- 2 £ ||y nl l(||y nl || < v n 6)\\ 2 < C6 



2-a 



D 



The next two results show that the centering constants in H4.6[) can be chosen in 
the same way as for (12. 2p when a^l. We say that a tempered stable law with index 
< a < 1 is centered if its log-characteristic function can be written in the form (12.51) 
where M is given by (13. ip . 

Theorem 4.4. If < a < 1 and v n satisfies (14.51) . then v~ l S n {n) =>- p x where p\ is 
centered tempered stable. 



Proof. In view of Theorem 14.31 and (I2.ip . it suffices to show that, if a n is defined by 
(14.71) . then a n — >• /riM^u xM(dx) . Fix < s < 1 and note that 



7? 



/ xP{v- x Y nX edx) + n / xP^Xi G dx) =: h + J : 

J{e<\\x\\<l} J{\\x\\<e} 
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Clearly, by (14. 9ft . lirrin^oo/i = Ji e< \\ a .\\ < nxM(dx) . Thus, it suffices to show 
(4.12) limlimsup ||/ 2 || = 0. 

Note ||/ 2 || < nv- l E[\\Y nl \\l(\\Y nl \\ < v n e)\ < nv- l [E(\\H\\l(\\H\\ < v n e)) + 
u n eP(||//'|| > v n e)\. Since a < 1, Karamata along with regular variation yields 
£(||#||l(||if|| < v n e)) ~ v n e 1 ' a P{\\H\\ > v n )a/{l - a). Then fl4TT2|) follows, using 
(14. 8 p and regular variation. □ 



Theorem 4.5. Suppose a > 1 and that for some (3 > a, 

(4.13) limsup sup u l ~^[a. — q(u, s)] < oo 

where v n satisfies (14. 5p . Then, t>~ 1 [S' n (n) — nE{H)\ =>- p 2 where p 2 is an infinitely 
divisible law with no Gaussian component, Levy measure M and mean 

(4.14) m=-a / j / (r - u)sU(du, s) 1 r- a-1 dnr(ds). 

Proof. Let 



oo 



(4.15) 6:= a / M (r-«)sn(d«,s) ^r" a_1 dr<r(ds) + / xM(dx). 

JO Js*- 1 {Jo J ./{||x||>l} 

We start with showing that the integrals on the right hand side of (14. 15ft are well 
defined. Let g(r, s) := / Q r (r — w)sll(dw, s) . It is easy to see that \\g(r, s)\\ < J^(r — 
u)H(du, s) =: g(r, s) . Clearly, 



f 
(4.16) g(j-,s)=r[l—ir(r,s)] + 

Jo 



d<K(u 7 s) 

u — du . 

ou 



and hence, 4pg(r,s) = 1 — 7i(r, s) = r a f°° u a 2+ ^u l @[a — q(u, s)]du = 
r a f^ u'~ a - 2+p u l ~^[a - q(u,s)}du + 0(r a ). Clearly (l4~T3j) holds with (3 replaced 
by (3 A 2. Thus, without loss of generality, we can assume that < 2. Define 
K = sup{u 1 ~^[a — q(u, s)] : s G S' ' -1 , < w < 1}. By hypothesis, K < oo. Thus, 

up- 1 y r 

where K' = Kj[a + 1 — /3) > since /3 < 2 and a > 1. Thus, as r 1 0, -§^g(r, s) = 
Oi^r 13 ' 1 ) uniformly in s, and hence for some C < oo, 

(4.17) g(r,s) <Cr p , r < 1, s G 5 d_1 . 
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It follows that f Q § sd _ 1 g{r ) s)r~ a ~ l dra(ds) < oo . It is easy to see that g(r,s) < r. 
Since a > 1, J x f sd -i g( r , s)r~ a ~ 1 dra(ds) < oo . Thus, the first integral in f |4.15|) 
is well defined. Since a > 1, it is easy to check that the second integral is 

ows, using Theorem 3.1.14 and Remark 3.1.15 in 



also well defined. Then it fol 



Meerschaert and Schefflerl (1200 ll ). that any tempered stable law with index a > 1 



has a finite mean. 

Next we want to show that 

(4.18) 



lim 

n— »oo 



/?. 



E(H) - a n 







Write nv- l E(H) - a n = nE \v~\Hr - Y nl )\ + nv~ l E [Y nl l(\\Y nl \\ > v n )} = h + h . 
Fix 1 < iV < oo and write 



n 



/ xP [y n 1 Y nl G dx) +n xP (v n x Y nl E dx) := I 2 \ + I22 ■ 

JiK\\x\\<N} Ju\x\\>N} 



'{l<\\x\\<N} 

By (J49J), it follows that 



(4.19) lim J 21 = / xM(dx) . 

n ^°° J{l<\\x\\<N} 

Using Karamata's Theorem we get \\I22W < nv~ l E [||F nl ||l(||F nl || > v n N)\ < 
nv^ElWHWldlHW > v n N)} ~ N l - a nP{\\H\\ > v n )a/(a - 1) as n ->• 00. This, in 
view of (14. 8p show that limAr^oo lim sup n ^ 00 H/22II = . In conjunction with (I4.19p . 
this shows that 



(4.20) lim I 2 = / xM(dx) 

J{\\x\\>l} 

It remains to show that 

/*oo f ( r r 

(4.21) limJ 1 =a/ / \ (r - u)sU(du,,i 
n ^°° Jo Js*- 1 {Jo 

To that end, fixO<e<l<iV<oo and note that 



-0-1 



dra(ds) 



h 



nE 



H 

W\ 



(«- 1 ||fT||-T)l(t;- 1 ||H||>T) 



n / / 9{ r i s )P \ v ii l \\H\\ ^ dr,—— ^ d s 
Jo Js d -i V II-" II 

+n [ g^ s )p( v -i\\H\\edr,-—eds 

Je Js d ~^ V 11-^11 



+n 



N 



/ TT 

g{r,s)P (v^WHW e dr,-r— 
-1 V \\-tt\ 



eds) =: In + In + Jig . 
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We shall now show that g is jointly continuous. Clearly, g(r, s) = g(r, s)s . Thus, it 
suffices to show that g is jointly continuous. Since q is assumed to be continuous in 
the second variable, an appeal to the dominated convergence theorem shows that tc is 
jointly continuous. By (jHJ), it follows that -§^^(u, •) is continuous for every u > 0. In 
view of (14.161) . it suffices to show that the function (r, s) i-4 J Q r u-^tt(u, s)du is jointly 
continuous. For that, fix a sequence r n —> r and s n — > s. Note that 

dn(u, s n ) dir(u, s) 



rn dir(u, s n ) , 
u du 



du 



r " dir(u,s) , 

u — - du + I u 

n ou J 



du 



du 



du =: J\ + J 2 



Clearly, as n — > oo, J\ — > f Q u-^ti(u, s)du . Let R = sup n>1 r n and note that, 



\Ji 



< 



R 



II 



dii(u, s n ) dir(u, s) 



R 



du du 

dir(u,s n ) dir(u,s) 



du< R 



d7r(u,s n ) dn(u,s) 



du 



du 



du 



du 



V 



du 



du + 2 



the second equality following from the identity \a — b\ = 2 (a V b) — (a + b). 
Since, |J^7r(w, s n ) V -J^tt{u, s)\ < —-^ir(u,s), an appeal to the dominated con- 
is continuous shows that 



da 



7T It, 



vergence theorem along with the fact that 

linin^oo J °° { J^7r(-u, s n ) V -^k{u, s) } du = — 1 , which in turn shows that J2 — > as 

n — > 00. This shows that g is jointly continuous. 

By ( 14. ip . (14. 8 p and the fact that g is jointly continuous, it follows that lim 



n— s-oo J 12 



-JV 



af f sd _ 1 g(r,s)r a 1 dra(ds) . Note that 



'■11 



< n 



/ / g(r,s)plv n 1 \\H\\edr,-—eds) 
Jo Js*- 1 V II-" II / 

< Cn [ [ r?P (v^WHW e dr, -^- G ds) 

Jo Js d -^ V ll-n|| / 



CnvJ 



r p P(\\H\\ edr)^C- 



G. 



-.0-a 



a(S' 



d-l- 



i (3 -a 

as n — > 00, using ( j4.17p . Karamata's Theorem, and ( 14. 8p . This shows that 
lim £ ^o lim sup n _ i>00 ||/n|| = 0. Finally, by similar calculations and the fact that 
Hfi 1 ^ s)|| < r i it follows that 

a 



/■oo 

\hz\\ <nv~ x \ rP(\\H\\ e dr) -4 

J Nv n 



-a{S d - l )N 



l-Q 



This shows that liniAr^oo lim sup n _ 
( 142TD . (I4TT8D follows. 



!13 



a — I 
0. Thus, I tOTjl follows. By ( IOOD and 
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From Theorem 14.31 and ( 14. 18ft it follows that 

77 77 

(4.22) v- l S n (n) - —E(H) = v?S n (n) -a n + a n - —E{H) =► p - 6 := p 2 



2001 



so that p2 has Levy representation [—9, 0, M\. Using (IMeerschaert and Schemer 
Remark 3.1.15), we can write the log-characteristic function of a tempered stable law 
with mean zero in the form (12. 4p . Then it follows easily that ( I4.14J) holds. □ 

Remark 4.6. As noted in Section 121 we can center to zero expectation in (I2.2p when 
a > 1, or dispense with the centering when a < 1. Theorems 14.41 and 14.51 shows that 
the same centering can be used for the tempered random walk. If a < 1, the limit is 
centered tempered stable, analogous to a centered stable law. If a > 1, and we center 
to zero expectation for the untempered random walk jumps, the limit contains a shift 
depending on the spectral measure and the tempering function. The shift comes from 
the fact that h = nv~ l E [H x - Y nl ] ->■ -m in ( fOTJ) . 

Remark 4.7. The speci al case d — 1 is also important in applications 



(IMeerschaert et al.l (120081 ) ). Suppose d = 1, and that H and vr(-, •) are as before. 
In this case, the conditional distribution of T given H can be written in a simpler 
form: P(T > u\H > 0) = k(u, 1) and P{T > u\H < 0) = n(u, -1). Let (v n ), (a n ) 
and p be as in Theorem 14.31 For n > 1, suppose that Y n i, . . . , Y nn are i.i.d. with 

Y nl = sgn(H)(\H\Av n T). 



Let S n (k) := J^ 7=1 ^nj- As a restatement of Theorem 14.31 we obtain that 

v^Snin) - a n ^ p. 
If a < 1, we can set a n = 0. If a > 1, we can take a n = nv~ l E(H), provided 



2a-q(u, 1) -<?(«, -1) 
hm sup 3— < oo 

ulO UP L 



for some (3 > a. 



Let {X(t)} be the Levy process generated by the tempered stable random vector 
X with distribution p, so that A(0) = almost surely, {X(t)} has stationary, inde- 
pendent increments, and X(l) = X in distribution. The next result shows that the 
tempered random walk (13.41) faithfully approximates the tempered stable process. 
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Theorem 4.8. Suppose that ( 14. 6 j) /ioWs as m Theorem \4.3\ Then 
(4.23) KS«([nt]) - io n } t > =► {X(t)} 4 > 

as n — > oo in the sense of finite dimensional distributions. 

Proof. The Levy representation of the limit p in (14. 6 j) is [0,0, M]. Use (14.91) 
to get [nt] P^- 1 ^ 6 •) ~ ntP{v~ l Y nl e -) -4- tM(-) and (OP]! to 
get lim^ohmsup,,^^ [nt]^ 2 ^||F n il(||F„i|| < w„5)|| 2 = 0. Then ^^([nt]) - 
Qn.t =>• /?/ follows by the ge neral convergence criteria for triangular arrays 



(IMeerschaert and Schemed . 1200 ll Theorem 3.2.2), where p t has Levy representation 
[0,0,£M], since 

a n t — [nt] / xP(v~ l Y n i G dx) < / ||a;||P(v^' 1 Yi l i G dx) 

•AnI<i) All*ll<i} 

< {v- 2 E\\Y nl l(\\Y nl \\<v n )\\ 2 } 1/2 ^0 



using (14.111) . To prove convergence of finite dimensional distributions, use the fact 
that increments of the random walk are independent. □ 

Remark 4.9. Take the exponential tempering function q(r,s) = e~ Xr for s = ±1. 
Then the random vectors X(t) have smooth density functions p(x,t) that solve 
a tempered fractional diffusion equation d t p = cqd°^ x p + c(l — q)d^' x p. where 
P(H < —r)/P(\H\ > r) ~ q as r — > oo. The operator on the right hand side 
is the negative generator of the continuous convolution semigroup associated with 
X. Some properties of th e temp ered fractional diffusion equation are developed in 



Baeumer and Meerschaert 



(120 10l ). Theorem 14.81 shows that the tempered random 
walk (13.41) provides a useful approximation to the process {X(t)}. In this case, the 
distribution of T, is given by P(Ti > u) = ir(u, s) = u a j r~ a ~ 1 e~ Xr dr, which in- 
volves the incomplete gamma function. The tempering thresholds v n do not depend 
on q. For example, if H belongs to the domain of normal attraction of some stable 
law, then we can take v n = cn l ^ a for some c > 0. Any random walk in the domain of 
attraction of a stable law can be modified using this tempering, to approximate an 
exponentially tempered stable. 

Remark 4.10. Suppose that the tempering variable is conditionally exponential with 
P(Ti > t\-Mn = s) := ir(t, s) = e~ Xst for some continuous s i— > X s > 0. Let h(r, s) = 
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r~ a ~ 1 q(r, s) and use ( 13. 2\\ to get u~ a e~ XsU = J h(r, s) dr. Take derivatives with 
respect to u on both sides to obtain -cra _a_1 e _As " — \ s u~ a e~ XsU = —h{u) and write 

(4.24) q{u, s) = u a+1 h{u) = {a + X s u)e- XsU . 

Using this tempering function for the Levy measure ( 13. ip yields a tempered stable 
law X with a particularly simple tempering variable Tj. If 1 < a < 2, then the 
form of the Levy measure shows that X is the sum of two independent exponentially 
tempered stable laws, one with index a, and the other with index a — 1. 

Remark 4.11. The goal of this paper is to construct random walk models that lead 
to a tempered stable limit. To conclude this paper, we provide a practical, heuris- 
tic interpretation of those results. A stable process serves to approximate a random 
walk with power-law jumps. A tempered stable approximates the same random walk, 
once the largest jumps are reduced. The tempering process represents an external 
force applied independently to each jump, the exact nature of which determines the 
tempered stable limit. Any random walk in the domain of attraction of a stable, and 
subjected to this type of independent tempering, can be fai thfully approximated by a 



tempered stable. A few concrete examples are provided in Meerschaert et al.l ( 120101 ): 



Precipitation data can be tempered due to atmospheric water content; measurements 
of hydraulic conductivity can be tempered by volume a veraging; dai l y sto ck returns 



could be tempered by automatic trading limits. See also lAban et al.l ( 120061 ) for addi- 
tional discussion. 
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